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Abstract 

The density of states of disordered systems in the Wigner-Dyson classes 
approaches some finite non-zero value at the mobility edge, whereas the 
density of states in systems of the chiral and Bogolubov-de Gennes classes 
shows a divergent or vanishing behavior in the band centre. Such types of 
behavior were classified as homogeneous and inhomogeneous fixed point 
ensembles within a real-space renormalization group approach. For the 
latter ensembles the scaling law p, = dv — 1 was derived for the power 
laws of the density of states p oc \E\^ and of the localization length 
£ oc \E\ ~ v . This prediction from 1976 is checked against explicit results 
obtained meanwhile. 

1 Introduction 

Some time ago I used real-space renormalization group arguments in analogy 
to the cell model of Kadanoff[T] in order to investigate the critical behavior [12] 
close to the mobility edge of the Anderson model[3]. Two types of ensembles 
were considered, a homogeneous and an inhomogeneous one. 

Homogeneous fixed 'point ensemble (HFPE). This ensemble is homogeneous 
in energy e. It is invariant under the transformation e — > e + constans. Since the 
density of states p stays constant during the renormalization group procedure 
the scale change 

r — > r/b implies e —5- eb d (1) 

with dimension d of the system. We assume one relevant perturbation to this 
system which grows like 

T^Tb V . (2) 

Depending on the sign of t the perturbation produces localized and extended 
states, resp. This perturbation is added to the HFPF in a strength increasing 
linearly in energy E 

t = cE, (3) 

where the mobility edge is taken at E — 0, and extended states at t > and 
localized ones for r < 0. c transforms under RG. 
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Inhomogeneous fixed point ensemble (IHFPE). In this ensemble the scale 
factors for length and energies are independent from each other. The ensemble 
is inhomogeneous in the energy, 

r ->■ r/b, e ->■ eb y . (4) 

It is assumed that there is no relevant perturbation to such an ensemble. 

Both ensembles yield power and homogeneity laws. The density of states 
obeys 

Phom = const., p inh oc fx = d/y - 1. (5) 

The localization length yields in both cases 

^k\E-E c \- v , v = l/y. (6) 

The low-temperature a.c. conductivity obeys the homogeneity relation 

, f b 2 - d a(ujb d ,Tby) HFPE . . 

^ T >-{b^a(^y,rby) IHFPE ' [) 

One deduces the d.c. conductivity in the region of extended states 

<7(0,7-)~7-», s={d-2)/y={d-2)v. (8) 

What comes out correctly on the basis of these ideas? Not only the scaling and 
homogeneity laws shown above can be deduced, but also such laws for averaged 
correlations, including the inverse participation ratio and long-range correlations 
between states energetically close to each other including those in the vicinity 
of the mobility edge. What has to be added are averages of matrix elements 
and of their powers for the transformation step by the linear scale factor b of 
the cell model. 

A short historical digression may be allowed. The oldest paper on the mo- 
bility edge i.e. the separation of localized and extended states of a disordered 
system was given by Phil Anderson[3] (1958) (Well aware of possible complica- 
tions by the Coulomb interaction he considered the transition from spin diffusion 
to localized spin excitations). It is a nice accident that its page number 1492 
coincides with the year of another important discovery. Earlier papers on dis- 
ordered systems, which became important for the development of this field was 
Wigner's[H |5] Gaussian matrix ensemble (1951) for nuclei. Probably the 
oldest paper on chiral systems is Dyson's paper[7] on disordered chains (1953). 
Other early contributions on disordered chains were by Schmidt |8 s and argu- 
ments that states in one dimension are localized. [9l 110). In 1962 Dyson gave 
the threefold classification of ensembles of orthogonal, unitary and symplectic 
symmetry depending on the behaviour under time-reversal invariance. [lll [T2"] . 

Since these early developments a lot of progress has been made. There 
are numerous calculations for the behaviour around the mobility edge both 
analytic and numerical. I refer to the review by Evers and Mirlin[13 . 1979 
marked important break-throughs: The scaling theories of localization by Abra- 
hams et al. [14] and by Oppermann and Wegner[T5] appeared. The mapping 
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onto a non- linear sigma-model was conjectured [IB], brought into its bosonic- 
replica 17 , its fcrmionic-replica [1 8] and finally in its supersymmetric[19j form. 
A self-consistent approximation for the Anderson transition was put forward by 
G6tze[30], Vollhardt and Woffle [21, 22]. A numerical renormalization scheme 
was devised by MacKinnon and Kramer [25]. 

Since then many more results and techniques were developped. Here I men- 
tion only a few: The complete classification of ten symmetry classes of random 
matrix theories, c-models, and Cartan's symmetric spaces was given by Zirn- 
bauer and Altland[23][25] and by Schnyder, Ryu, Furusaki, and Ludwig[5B] after 
several occurences of chiral and Bogolubov-de Gennes classes [271 1251 1291 |3"01 I3T1 
[521 13"51 IM] . These classes are listed in table Q] since I will refer later to this 
nomenclature. 

Table 1: Symmetry classes of single particle Hamiltonians defined in terms of 
presence or absence of time-reversal symmetry (TRS) and particle-hole symme- 
try (PHS). Absence is denoted by 0, presence by the symmetry square ±1. SLS 
indicates absence (0) and presence (1) of sublattice or chiral symmetry. After 
ref.126] 



System 


Cart an 


symmetry 


TRS 


PHS 


SLS 




nomenclature 










standard 


A 


unitary 











(Wigner-Dyson) 


AI 


orthogonal 


+1 










All 


symplectic 


-1 








chiral 


AIII 


unitary 








1 


(sublattice) 


BDI 


orthogonal 


+1 


+1 


1 




CII 


symplectic 


-1 


-1 


1 


Bogolubov- 


D 







+1 





de Gennes 


C 







-1 







Dili 




-1 


+1 


1 




CI 




+1 


-1 


1 



Transfer matrix approaches originally used for linear chains were developped 
for the non-linear q-model[19j as well as for the for the distribution func- 
tion of the transfer matrix of chains with many channels (DMPK-equation[35, 
36J). These techniques allowed the determination of correlations, wave- function 
statistics and transport properties. Such chains can have broad distributions of 
conductivities and even cases of perfect transmissions were found [37] 138] . 

In two dimensions some of these symmetry classes allow the inclusion of a 
topological 0-term. As observed by Pruisken et al. (39] [40] the Wigner-Dyson 
unitary class with this term describes the integer quantum Hall effect. Another 
term which may be added is a Wess-Zumino term. Such terms are of impor- 
tance in the study of disordered Dirac fermions, which appear in dirty d-wave 
superconductors [1H H2j [43] and in disordered graphene. |S1 H3 146 ] |47 [ 148 ] 149] 

Network models originally introduced by Shapiro [50] are very useful for the 
description of quantum Hall systems as first shown for the integer quantum Hall 
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effect in the Chalker-Coddington-model[5"T] . 

Obviously the HFPE applies to the Wigner-Dyson classes, whereas the 1HFPE 
applies to the chiral and the Bogoliubov-de Gennes classes. 

The main object of this paper is the comparison of the scaling law for the 
IHFPE 

H = dv-l. (9) 

derived from ([SJ [5]) with various results meanwhile obtained. I will shortly 
comment on the scaling law ([8J for the conductivity in subsection l4.ll 



2 One dimensional chains 
2.1 Thouless relation 

Thouless[52] following Herbert and Jones [53] considered a one-dimensional chain 
governed by the Hamiltonian 

N N-l 

H = 5^ei|»><i| - ]T (Vi, i+1 \i)(i + 1\ + V i+hi \i + (10) 

i=l i=l 

and found in the limit N — > oo that the function 

K(z) — J dxp(x) hi(x — z) — In |V|, — 7r < arg ln(a: — z) < it (11) 

connects both the integrated density of states 1(E) and the exponential decrease 
of cigcnf unctions X(E) (inverse correlation length £) 

K(E±iO) = X(E) Tinl(E). (12) 

The density of states is symmetric in chiral and Bogoliubov-de Gennes classes 
p(-E) = p(E). Then besides K(z*) = K*(z) also K(-z) = K(z)+ms(z) with 
s(z) — signQ ! (z) holds. If K (z) + ms/2 oc z 1 for small z, then 

K(z) + i y = c ^e i7(<# - S7r/2) , z = re i<f> (13) 

with real c. Then 

K(E + iO) = c\Ep (cos(| 7 ) - i sign(£;) sin(| 7 )) - |i, (14) 

from which A oc |-E| 7 , p(E) oc |i?| 7_1 follows in agreement with ([9]). One 
observes that 

^=/d^. (15) 
az J z — x 

Thus the sign of the imaginary part of (ITBl is opposite to the sign of 3z. This 
implies that 7 < 1. If a contribution with 7 > 1 appears, then there is also 
a contribution with 7=1, which according to (| 14[) does not contribute to A, 
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but to a finite density of states in the center of the spectrum. Such a system is 
described by the homogeneous fixed point ensemble. 

For 7 < the integrated density of states would diverge. Thus these argu- 
ments can only be applied for < 7 < 1. 

In a number of cases the asymptotic behaviour is given by a power multiplied 
by some power of the logarithm. Then 

K{Z) + iy = cr 7 e i7W>-W2) ( lnr + i(0 _ s7r /2))9. (16) 

This yields for 7 = and 7 = 1 

_j K(E + iO) 

~0 c(ln \E\)a - icgf sign(£)(ln l-EQf- 1 (17) 
1 -icE(ln \E\)3 -cgf I^Klnl^Df^ 1 

7 A ~ p ~ 



and thus 



(ln\E\)e ^f^- (18) 

1 I^Klnl^Df- 1 (ln\E\)s 

Dyson[7] calculated the averaged density of states for the chain (|10p with 
e.i = and random independently distributed matrix elements V, for which he 
assumed a certain distribution and obtained 

^ ~ MinW (19) 

which corresponds to the case (fT8")l with 7 = and g = — 1. Indeed Theodorou 
and Cohen[S3] and Eggarter and Ridinger[SS] found the averaged localization 
length diverging 

£HH£|I (20) 

in agreement with ([TJ 



2.2 Ziman' s model 

Ziman 56 (compare also Alexander et al.|57j) considered a one-dimensional 
tight-binding model (his case II) ([TO]) requiring the diagonal matrix elements to 
vanish, ej = 0, and the hopping matrix-elements to agree pairwise V2 m ,2m+i — 
V2rn+i,2rn+2- Apart from this restriction he assumed the Vs to be independently 
distributed with probability distribution 

p (V) = (1 - a)V~ a , 0<V<1, -oo<a<l. (21) 

He obtained for these distributions 



-1< a < 1 



v 

2(l-q) 



3— a 3 — o 

-3<a<-l i=a 
a < -3 2 



(22) 



Obviously the first row describes models in accordance with the IHFPE, the 
second and third row with the HFPE. 
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2.3 Further one-dimensional results 

Titov et al. [5S| have summarized and completed results for the density of states 
of all classes of chains with N channels as shown in table [5J 

Table 2: Density of states close to E — for various universality classes (after 



Titov et al.[55]) 


Class 




p(£7) x = .Br 


Chiral 


all classes, odd N 


Ix-^n-^x)) 




kill, even N 


|xlnx| 




CII, even N 


|x lnx| 




BDI, even TV 


lnx| 


BdG 


CI 


W 




C 


x 2 




D,DIII, N + 2 


\x~ x lrr 3 {x)\ 




D,DIII, N = 2 


two mean free paths 



All chiral classes are equivalent by a gauge transformation for N = 1 and 
yield the Dyson result (fT§j) and £ oc | In \E\ | for this case in agreement with (ITS)) . 
Due to Gruzberg et al.[S2] also the BdG classes BD and Dili fall into the same 
universality class. The localization length docs not diverge for the chiral classes, 
if N is even. The same holds for (general N) for the BdG classes C and CI. 

3 Bosons From One To Two Dimensions 
3.1 One-dimensional chain 

Whereas the Hamiltonian (|10l) yields the equation for eigenstates = 

Eipi = eiipi - V iti -iipi-i - V i>i+ iipi + i (23) 

one obtains a similar equation for harmonic phonons governed by the Hamilto- 
nian 

i i 

which reads 

mu 2 Xi = Wi(xi - x i+ i) + Wi-!(xi - Xi_i). (25) 

Thus Thouless' arguments can be applied again with x = to 2 . Since there are 
no states for lu 2 < 0, one has 

K(z) = crT e i7W-«r) j z = re i ' t> , (26) 

which yields 

X(w 2 + 10) = cuj^e' 1 ^ (27) 



G 



and thus 

A(w) = cw 27 cos(tt 7 ), I(uj) = cw 27 sin(7T7), < 7 < 1. (28) 

Alexander et al.[57] (cases a,b, and c) and Ziman[55] (case II) determined the 
density of states p(u>) oc uj^ for harmonically coupled phonons with indepen- 
dently distributed spring constants 

p(W) = (1 - a)W- a . (29) 

Ziman moreover determined the localization length £ oc lo~ v and obtained 





V 




< a < 1 


2(1-") 


a 


2-a 


2-a 


-1 < a < 


1 -a 





a < -1 


2 






(30) 



Again the first line is in agreement with the IHFPE, whereas the two other cases 
correspond to the HFPE. 

3.2 Bosonic excitations discussed by Gurarie and Chalker 

Gurarie and Chalker [60 point out that bosonic systems with and without Gold- 
stone modes show different localization behavior. 

John et al. [61] investigated localization in an elastic medium with randomly 
varying masses. For d > 2 they found extended states for small frequencies 
uj. The phonon states are localized beyond some critical uj c . This transition is 
described by the orthogonal ensemble. For d < 2 all states are localized and 
obey [61] v = 2/(2 — d). The density of states for phonons shows the same power 
law p(uj) oc as in the ordered case. In this system with Goldstone modes 

the critical density below which the density of states would differ from that of 
the ordered system is d c = 0. 

In a disordered antiferromagnet one obtains below the critical dimension 
d c = 2 

pM-uA H=^£, ««)~^ ( 31 ) 

where the result of [62] and the argument of[60 have been generalized from 
d = 1 to general d < d c . This is in agreement with the IHFPE. These results 
rest on the assumption that there is a single length scale £ oc 

4 Electronic Systems In Two Dimensions 

4.1 Conductivity in two dimensions 

From the homogeneity law s — (d — 2)u, ©, which works well for d > 2, I 
concluded [2^ s = for dimensionality d = 2 and thus a jump to a minimum 
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metallic conductivity. At that time I did not expect that v may diverge as d ap- 
proaches 2. This was found three years later by means of explicit renormalization 
group calculations [31 115) . Thus the idea of a minimum metallic conductivity 
was in error for the orthogonal and unitary Wigner-Dyson classes, where all 
states are localized in d = 2. The critical conductivity in the symplectic class 
shows some distribution [55] and is of order e 2 /h. 

In many two-dimensional models of chiral and Bogolubov-de Gennes classes 
including the classes applying to d-wave superconductors and graphene one 
obtains a finite non-zero conductivity of order e 2 /h at criticality. This is to a 
large extend due to edge currents, as first observed by Pruisken et al. [SHIED] f° r 
the integer quantum Hall effect. Thus although the prediction^ turns out to 
be correct, the true mechanism is more complex. 

4.2 Chiral and Bogolubov-de Gennes models in d = 2 di- 
mensions 

The unitary case of chiral models (Gade and Wegner, Gade[3_U [32]) yields at 
intermediate energies effective exponents 

12 

y = -g, /x = -i + -g (32) 

in accordance with At asymptotically low energies p oc E^ 1 ^ 2 (E) corre- 
sponds to the limit B — > oc. These results as well as many similar results for 
various disordered Dirac hamiltonians are obtained under the assumption that 
the localization length is given by the cross-over length from chiral to Wigner- 
Dyson behaviour without taking further renormalization into account [3"2"j ; see 
also the argument after eq. (6.60) of the review by Evers et MirlinJUJ. Alter- 
natively the integrated density of states from the band center up to energy E 
is set to £~ 2 as in Motrunich et al.[51], which yields (JSJ per definition. It is 
important that only one coupling yields a relevant perturbation. The conduc- 
tivity itself stays constant for the chiral models in d — 2. The exponent which 
drives the renormalization of the energy is usually called dynamical exponent z, 
which is identical to the exponent y of ([9]) . A more rigorous investigation of the 
localization of such systems taking into account any dependence of the initial 
couplings on the energy and of the cross-over would be of interest. 

The spin quantum Hall effect yields at the percolation transition point [65l 

u = 4/7, M = l/7 (33) 

in agreement with ([9]). The same behavior is obtained for the Bogolubov-de 
Gennes class C if two of the four nodes of a dirty d-wave superconductor are 
coupledgUiS]. 

4.3 Power law for density of states, finite localization length 

The two fixed point ensembles describe the situation, in which the localization 
length diverges and the density of states either approaches some finite non- 
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zero value (HFPE) or diverges or goes to zero by a power law, which may be 
augmented by a logarithmic term. As mentioned above this holds for chains 
with an even number of channels in the chiral classes and for the Bogolubov-de 
Gennes classes C and CI. Certain single-channel models of class D and Dili 
show a divergence of the density of states p oc \E\~ 1+S without divergence of 
the localization length [55], 

Gurarie and Chalker|69) found that bosonic excitations in random media, 
which are not Goldstone modes, obey paw 4 with finite localization length at 
low frequencies. 

Apparently this type of behavior is not covered by HFPE and IHFPE. 

5 Conclusion 

The scaling prediction © of the IHFPE relating the exponent of the density of 
states and of the localization length yields correct results in the cases, in which 
I found both exponents. The author appreciates the wealth of systems, which 
has been found and investigated over the years. 
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